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>■ ' The strong (yf-coupling characterizes the interaction of heavy mesons with 

00 ' 



pions in typical vertices H*H7i, H*H*7r, where [H*; H) stands for vector and 
pseudoscalar {B*;B) or {D*;D) heavy mesons. Its estimation by different 
theoretical methods has led to a wide range of possible values. We describe 



p i' a new approach to the determination of g, which exploits the rare radiative 

decays B* ^ B'-fy and D* ^ D'-f-f. It is shown that the branching ratio of 



/\ ' D* -^ D'j'j can be expressed as a function of a single unknown g and we 



calculate it to be in the measurable range between 1.6 x 10 ^ and 3.3 x 10 ^ 

for 0.25 < ^ < 1. 

PACS numbers: 12:39 Fe; 12.39 Hg; 13.25.-k; 13.40 Hg 
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1. Introduction and Experimental Overview 

We consider here the strong and electromagnetic decays of the heavy 
vector mesons D* and B* of spin-parity 1^, with the special aim of gaining 
information on their strong couplings gB*Bn, gD*Dw It is well known that 
the main decays of D*'s proceed either as a strong transition D* — ^ Dir 
with a final pion of about 40 MeV/c momentum, or as an electromagnetic 
transition D* — > D'j with a photon of momentum of nearly 140 MeV/c. On 
the other hand, since the mass difference Mb* — Mb is only 45.8 MeV, the 
decay B* -^ Bn cannot take place and the main decay of B* is the radiative 
process B* -^ B'j. 

The interaction between mesons containing a single heavy quark Q and 
pseudoscalar Goldstone bosons is presently best described by an effective the- 
ory [1,2,3] which contains flavour and spin symmetries in the heavy mesons 
sector and chiral SU(3)l ® SU{3)ji symmetry in the light one (for a recent 
review, see [4]). As a result of the heavy quark symmetry, the heavy me- 
son chiral Lagrangian {HMxL) which implements this scheme contains just 
one coupling, denoted by g, to characterize the strength of several vertices, 
0*0*71, 0*071, B*B*7i and B*B7v. While the gB*BTT vertex is not accessible 
in a direct decay, the gD*Dw coupling may be measured in principle from the 



decay width of the decay channel D* -^ Dn. So far, only an upper limit 
is available from the ACCMOR Collaboration at CERN [5], T{D*+) < 131 
KeV based on the high-resolution measurement of 127 D*^ events. 

In view of the great interest in the strength of the B*B7i and D*D7r 
vertices, which are relevant for the analysis of various B and D decays, a 
large number of calculations has been performed to obtain the (^-coupling. 
The results are quite divergent, as it will be described in the next chapter. 
Recently, Daphne Guetta and myself [6] have suggested a new approach to 
the determination oi g. We considered the decays B* -^ B'-f-f and D* -^ D'-fj 
within the framework of HMxL and we have shown that the measurement 
of the branching ratios and spectra of these decays can provide information 
on the (yf-coupling. This method is especially profitable in the charm sector. 

As a background for the model we have built, I review here succinctly the 
experimental status of the B* and D* decays and present a summary of the 
theoretical attempts to describe the main decays to a pion or a photon. 

The main decay of B*, the electromagnetic transition B* -^ B'-f has been 
observed both at the Cornell Electron Storage Ring (CESR) [7] and at LEP 
[8]. There is obviously no measurement of the width of this transition. It has 
been studied in a variety of theoretical models including quark models [9], 



the chiral bag model [10] followed by effective chiral Lagrangian approaches 
[11], potential models [12], QCD sum rules [13] and an analysis of experi- 
mental D* branching ratios using HMxL [14]. The predictions range from 
T{B*"{B*+) -^ B%B+)-f) = 0.04(0.10) KeV [13] to 0.28(0.84) KeV [10,11], 
with most calculations concentrating in the higher range. 

The D* meson was discovered more than twenty years ago [15] and its 
decay branching ratios have been studied extensively. The current PDG av- 
erages [16] are Br(D*+ -^ D+n") : Br(D*+ -^ D°ti+) : Br(D*+ -^ D+7) = 
(30.6 ± 2.5)% : (68.3 ± 1.4)% : (1.1 ± 2.1 ± 0.7)% and 5r(D*° -^ D^n") : 
Br(D*° -^ D^-f) = (61.9 ± 2.9)% : (38.1 ± 2.9)%. A recent [17] CLEO ex- 
periment on D*^ decays gives the more accurate branching ratios By{D*^ — > 
D+n") : Br(D*+ -^ D"tt+) : Bt{D*+ -^ D+7) = (30.7 ± 0.7)% : (67.6 ± 
0.9)%: (1.7 ±0.6)%. 

These D* electromagnetic and strong decays have been calculated with 
the same models used for B* and most papers of Refs. [9]- [14] have considered 
also D* decays. Additional calculations referring to D* only include quark 
models [18], the chiral bag model [19] and use of sum rules [20]. Here again 
the theoretical calculations span an order of magnitude range for the predic- 
tion of the widths, from r(D*°) ~ (3 - 10)KeV [13] to (60-120)KeV [9,11]. 



Many of these theoretical results are fitted to obtain correct relative widths. 
The real test will come when the absolute width will be measured. Accord- 
ing to several of the models, typical values are r{D*^ — ^ all) ^ 80KeV, 
r[D*° -^ all) ~ 60 KeV [9,18,19], not far from the present upper limit [5]. 



2. The Heavy Meson Chiral Lagrangian and Estimation of g 

The treatment of the physical processes involving soft pions we described 
in the previous section is performed within the framework of an effective 
theory, the "Heavy Meson Chiral Lagrangian" (HMxL) which embodies 
two principal symmetries of Quantum Chromodynamics. At one end, there is 
the SU(2Nj) heavy flavour-spin symmetry characteristic of the infinite heavy 
quark mass limit. In this limit, the interactions with the light hadrons are 
independent of the mass of the heavy quarks; moreover, the independence of 
the interaction on the spin sq of the heavy quark allows to define degenerate 
doublets of heavy states with spin-parity Sq = {se ± |)^ where sg is the spin 
of the light quark, which in our case is g = -u, d, s. 

The doublet members are the pseudoscalar and vector mesons correspond- 
ing to s^ = I and we make the assumption that both the c and the b quark are 
sufficiently heavy. Symmetry-breaking corrections to this scheme are taken 
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into account by terms obtained in a -p— expansion [4]. 

At the other end of the energy scale is the chiral hmit of QCD real- 
ized for M^ -^ 0(£ = u,d,s); the QCD Lagrangian is then invariant under 
SU(3)l X SU(3)/j transformations. This symmetry is spontaneously broken 
to the vector subgroup SU(3)y and the resulting Goldstone bosons are the 
eight pseudoscalar mesons n, K and rj. The quark mass terms breaking the 
chiral symmetry help the Goldstone bosons to acquire their mass. 

The effective Lagrangian which contains these symmetries is expressed 
[1-4] in terms of the hadronic fields for heavy and light mesons. The heavy 
vector {B*,D*) and pseudoscalar {B,D) mesons are represented by a 4 x 4 
Dirac matrix H, with one spinor index for the heavy quark and the second 
one for the light degree of freedom, 
1 + 



H 



PX -Pi5 , H = ^Mho . (1) 



P* and P are the respective annihilation operators of vector (1^) and pseu- 
doscalar (0^) heavy mesons with four- velocity v^. 

The Goldstone bosons are represented with the aid of a unitary 3x3 
matrix X) = exp(2iM//) with the M being the usual 3x3 hermitian traceless 
matrix describing the octet of pseudoscalar bosons. 

The most general Lagrangian describing the interaction of heavy mesons 
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with Nambu-Goldstone bosons, which is invariant under Lorentz transfor- 
mations, parity, heavy-quark spin flavour symmetry and chiral symmetry is 
given in the framework we described by [1-4] 

C = iTr {SaV^D^baHt} + ^d^J^abd^^l 

tgTr{H,-f,-f,A':,Ht,} (2) 

where D^ is a covariant derivative. D^ = d^ + V^, and / is the pion decay 
constant, / = 132 MeV. 

Using E(x) = ^^{x) one introduces a vector V^ = |('C^f^M^~^^M'^^) ^^^ ^^ 
axial current A^ = ^{^"^d^^+^d^^"^). The first two terms of (2) are the kinetic 
energy terms and the third one is an interaction term, defining the strong- 
interaction coupling g. Explicit interaction terms are obtained by expanding 
the axial current in (2) and keeping the first term A'^ = {i/ f)d^M + . . .. 
Expressing the interaction in terms of {D*,D) fields (the same holds for 
B*,B fields) one has 



£int 



eflf 



—Dld^'MD^ + h.c 
f ^ 



+ ^e^,.,Z}*^9'^MD*t-^r (3) 



and we see that, e.g., D*D7i and D*D*7i vertices are characterized with the 
same strength g. The (^-coupling is directly related to the hadronic coupling 



gD*Dn which is defined by the on-shell matrix element 

{D"{p)7r+{q)\D*+{p + q)) = gn*D.e,q>^ (4) 

where e'^ is the polarization vector of D*^. Likewise, 

(D*°(p,ei)7r+(g)|D*+(p + g,e2)) = ^D-D^.e^-re^aPV • (5) 

/^From (3-5) one finds the link 

gD'DTT = QD'D'n = —^9 ■ (6) 

Moreover, isospin symmetry requires 

Qd'Ott = gD*+D°iv+ = — V 2(7£i.+£i+7ro = y2g£)'"'D°TT° = —9d*°d+tt- ■ (7) 

The interest in the value of g is not limited to the theoretical interest 
in the strength of the axial interaction defined in Eq. (2). The knowledge 
of g is of great phenomenological value, since its strength is required in the 
analyses of many electroweak processes [4]. Among these, for example, are 
heavy-to-light semileptonic exclusive decays like B -^ Triu, Dg -^ Kiv which 
are promising processes for the extraction of CKM matrix elements like 
iKtfel, and their analysis with VMD requires the knowledge of gB*B-K- Other 
processes requiring such knowledge are B -^ D*Tt(!.v decays, chiral corrections 



to B ^ D processes, decay constants of heavy mesons, radiative processes 
like D* -^ D^, B* -^ B'-f and more. Thus, no wonder that a large number of 
theoretical papers has been devoted during the last years to the computation 
oi g. In the rest of this Section we present a succint overview of the main 
attempts in this direction. 

The theoretical attempts may be grouped into several classes; however, 
even within a single class of models, the variation of g turns out to be quite 
large. The non-relativistic quark model leads to [3] the largest value g = 1, 
while slightly modified quark models [18, 21] bring this value down to about 
g = 0.8. On the other hand, in a calculation [22] in which the effect of 
the relativistic motion of the light antiquark is taken into account by the 
use of Salpeter equation one arrives at g = |. Somewhat higher values 
were obtained in recent quark-model calculations: with a relativistic quark 
model based on the light front formalism Jaus obtains [9] g = 0.56; in a rel- 
ativisitic quark model with direct quark-meson interactions [23] one arrives 
ai g = 0.46; and a quark model with Dirac equation [24] finds g = 0.61. The 
QCD sum rules have also been used [25] extensively to the calculation of 
g. The outcome of these calculations is generally in the direction of small g 
values, between 0.15 and 0.35, which would imply that the decay width of D* 
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is rather small, i.e. below 45 KeV. Finally, we mention a recent lattice QCD 
determination [26] oi g = 042(4) (8) and the analysis of Stewart [14] of the 
experimental data on D* -^ Dn, D'-f which incorporates symmetry breaking 
terms in the Lagrangian and deduces g = 0.27lo;o4. 



3. A Model for Two-photon Decays D* -^ -D77, B* -^ B-yj 

Recently we have considered the two-photon processes B*{D*) -^ B{D)'y'j 
[6], not discussed previously in the literature, by using the HMxL and we 
found that the D* -^ D77 could provide a measurement of the much sought 
after (yf-coupling. 

The calculation of radiative processes requires the addition of the electo- 
magnetic interaction to the Lagrangian of Eq. (2). This is performed by the 
usual procedure of minimal coupling which leads to the replacement of deriva- 
tive operators by covariant derivatives containing the photon field. However, 
this does not suffice to account for the observed magnetic dipole transitions 
B* -^ By, D* -^ Dy; to account for these, a contact gauge invariant elec- 
tromagnetic term proportional to F^j^ must be added, which has the form in 
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the heavy mass hmit [11, 14] 



£(/^) = ^Tr{H,a,,F^''H,6at), (8) 



where /x is the strength of this anomalous magnetic dipole interaction and 
has mass dimension [1/M]. 

In the following, we present our calculation for the charm sector; then, 
we shall comment on the features arising in the beauty sector. From (8), 
additional electromagnetic vertices obtain representing D*D*'j and D*D'y 
interactions of same strength. These are 

{^{k,e)D{vi)\D*{v2,e2)) = -teMn^fie^^^fse'^k'^v^el (9) 

(7(A;, e)D*{vi, ei)\D*{v2, ea)) = e/iMz)*(ei ■ ke ■ €2 - €2 ■ ke ■ ei) . (10) 

The calculation is performed to leading order in chiral perturbation the- 
ory and to this order there are no counterterms [14,27]. In addition to the 
Feynman diagrams obtained from Eq. (2) with minimal electromagnetic in- 
teraction and from Eq. (8), we must include to the same order the pion axial 
anomaly, whose strength is known. All these terms are of the same order in 
an 1/Nc expansion. 

We start with the description of the decay of the neutral D*. The decay 
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amplitude A for D*° -^ D°^^ may be written as 

6 
A = ^anomaly + ^tree + /_^ ^loops • (-'--'-J 



1=1 



We shall describe now the eight contributions to the amplitude, indicating 
the couplings entering into each of them, without giving here the detailed 
expressions which can be found in Ref. [6]. 

^anomaly represents D*° -^ D°"-k" -^ D^'yy via a virtual neutral pion. 
Since the physical decay D*° -^ D°tt° is allowed, we limit ourselves to a 
region for s = (/ci + A;2)^ which goes up to 20 MeV away from the pion 
mass. Given the strength of the pion axial anomaly of (a/vr/), where / is 
the pion decay constant and a = e^/An, Aanomaiy is proportional to ago'Dn- 
The tree level graph is due to the transition D*° -^ iijj*o n^ _^ -D°77, 
containing two insertions of the anomalous magnetic operator. Hence, A^^^e is 
proportional to a/i^. A\J describes the transition D*" -^ (D*^n^) -^ D°77 
with the two photons radiated from the virtual charged pion. Additional 
graphs, required by gauge invariance have one photon radiated from the loop 
and the second emitted from the D*D*7i, D*D7t vertices or both photons 
emitted from these vertices. The other loop diagrams, ^loop ^ ^loop come 
from diagrams where both the strong coupling and the magnetic one are 
involved. A^^^ is given by D*° -^ "L)*°"7 -^ {D*+n-)-j -^ 1^°77, being 
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(3) 

thus proportional to agD*D*TTgD*DTrl^- In ^loop ^ D*"D°'y vertex replaces the 
D*°D*°'~f one in the initial step, the diagram being proportional to agj^tD^fi. 
A\J^^ describes the transition D*° -^ (D^n^) -^ D*"^ -^ -D°77, where the 
first photon is emitted by the charged pion from the (D+vr^) loop, while the 
second one comes from the D*° — * D°'~f transition. The expression is then 
proportional to ag'j^^j^^ii. Exchanging D with a D* in the loop one gets A\J^^, 
which is thus proportional to ag'jj*j^,^^. Finally, we have a diagram given by 
D*° -^ (£)+7r^) and then the virtual pion emits one photon while the virtual 
D^ emits the other one becoming D*^ . The virtual (D*+7r^) recombine to 
a D". This is described by A\J^^ and is proportional to agjyt£)^jj,''^\ where 
/i'^+) is the strength of the D*^ -^ D+7 transition. 

Calculating from (11) the decay width, one obtains [6] an expression 
containing 13 terms, which depend on various products g" fi^ fi^^''' ^ where 
a,j3 have values to 4 and 7 has values 0-2. At this point, the crucial step 
is to use the existing experimental information on the relative branching 
ratios V{D*° -^ D^n") : T{D*° -^ D°j) = (61.9 ± 2.9)% : (38.1 ± 2.9)% and 
T{D*+ -^ £)%+) : T{D*+ -^ D+7) = (67.6 ± 0.9)% : (1.7 ± 0.6)% [16,17]. 
This allows us to establish /i ^ Q.Qg/Mo* and /i+ ~ l.Tg/Mo*. As a result. 
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we were able to express r{D*° -^ D°'~f'~f) [6] as a function of g only: 

T{D*" -^ L)°77) = [2.52 x 10"^^^^ + 5.66 x 10"^^/ 

+ 4.76 X lO^y + 3.64 x 10"^°^^ + 1.53 x lO^V] GeV . (12) 

In obtaining (12) we assumed that the couphng constants are relatively posi- 
tive, as indicated by theoretical analysis [14]. However, if we assume opposite 
sign for various pairs of couplings, we found that the changes are rather small, 
the reason being that the main contribution is given by quadratic terms. 
Let us define now the branching ratio for this decay 

Bfl(B-"^B^.) = ^^j^;^^l,^|,jL__. (13) 

In the denominator, we can use for r{D*° -^ D°n°) the expression obtainable 
from Eq. (4) 

r(Z}- ^ DV) = -^f^lp.l' = 1.25 X lO-y GeV (14) 

while for T{D*" -^ D^-f) we use the experimental fact [16] T{D*" -^ D"n") : 
T{D*° -^ D°-i) = 61.9 : 38.1, to rexpress T{D*° -^ D°-f) in terms of g^ as 
well. Thus one obtains [6] 

T{D° -^ all) = (2.02 ± 0.12) x 10"^^^ GeV (15) 
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and as a consequence 

^ ^ ^^^ 2.02 X 10-4c/2 ^ ^ 

Obviously, a measurement of this ratio will constitute a measurement of g. 

Turning to the B*" -^ B"'~f^ decay, one has a rather different situation. 
Firstly, the branching ratio is now defined as 

BRiB*'' -^ B°-i-i) = -^ ^-^ (17) 

since there is no strong decay of B*°. The quantity in Eq. (f 7) turns out to 
be a function oi g, fi and fi~^ and an analysis similar to the one we performed 
for D*° decay is not possible. Still, a certain amount of information is ob- 
tainable by a judicious analysis in the above parameter space [6]. However, 
we shall not address this topic here. 



4. Discussion and Summary 

Firstly, some remarks about the framework of the calculation of [6] . The 
analysis of D*° -^ D°'-yj has been performed to leading order in chiral per- 
turbation theory and mostly to leading order in a 1/M expansion. Correc- 
tions to the leading order of (2) have been studied extensively in recent years 
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[4,14,28]. A comprehensive treatment of such corrections is beyond the scope 
of the calculation presented in Ref. [6]. However, one notes that several fea- 
tures belonging to the next order are included in their [6] treatment, like 
the use of physical masses for the degenerate doublet of heavy mesons in the 
decay calculations and in propagators. For the latter, as well as for vertices 
and normalizations the convention of Ref. [4] is used. Thus, the propagator 
of heavy vector mesons is given by —i{g^" — v^v'^)/2[{v ■ k) — A/4] and of 
pseudoscalar mesons by i/2[{v.k) + 3A/4], where A = Md*{Mb*) — Md{Mb) 
and V, k are the velocity and the residual momentum. 

Additional technical points, which should be mentioned are: (i) the loop 
calculations for D*° -^ D°'~f'-f include also contributions from intermediate 
states containing i^— mesons, like D*° -^ (K^D*^) -^ D°'-yy with the pho- 
tons emitted from the virtual K~'s; (ii) contributions from diagrams contain- 
ing three heavy meson propagators were neglected, as these are very small 
indeed; diagrams with two heavy meson propagators were included; however 
their contribution is quite small; (iii) the contribution of the ?7°-anomaly 
has been estimated and found to be small; (iv) the ofF-the-mass-shell q^- 
dependence of the anomaly has been neglected. 

The calculation described here was performed for neutral D*° decay (like- 
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wise for B*). Obviously there are also the D*^ -^ D+77, B*^ -^ B^'-fj 
and D*^ -^ -D+77 decays. For these decays, one has to consider also the 
bremsstrahlung radiation emitted by the initial or final charged particles. 
We have estimated these decays and we found that the bremsstrahlung part 
is orders of magnitude larger than the direct one; hence a different type of 
analysis is required [29] and we do not address this here. 

To summarize, we have shown [6] that the measurement of the D*° — > 
D°'~f'~f branching ratio constitutes an ideal tool for obtaining the magnitude 
of the strong (jf-coupling, since it can be expressed in terms of g only (Eq. 16). 
This has been achieved by combining a theoretical calculation of D*° -^ D°'y 
using the Heavy Meson Chiral Lagrangian with the experimental information 
relating the electromagnetic {D* -^ D'-f) and strong {D* -^ Dn) partial decay 
channels. The various theoretical estimates put g in the range 0.25 < g < 1. 
Hence, from Eq. (16) one obtains 

g Br[{D*" -^ D^-f-f) / D*" -^ all] 



0.25 


1.7 X 10-6 


0.38 


3.9 X 10-6 


0.5 


6.9 X 10-6 


0.7 


1.4 X 10"^ 


1 


3.3 X 10-5 



These figures indicate that the suggested measurement is indeed feasible 
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in the not too distant future. 
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